On the consistency of warm inflation 
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Conditions are obtained for the existence of a warm inflationary attractor in the system of equa- 
tions describing an inflaton coupled to radiation. These conditions restrict the temperature depen- 
dence of the dissipative terms and the size of thermal corrections to the inflaton potential, as well 
as the gradient of the inflaton potential. When these conditions are met, the evolution approaches a 
slow-roll limit and only curvature fluctuations survive on super-horizon scales. Formulae are given 
for the spectral indices of the density perturbations and the tensor/scalar density perturbation 
amplitude ratio in warm inflation. 



I. INTRODUCTION 



o , 

Inflationary models 0, H, 0] have proved very sucessful in explaining many of the large scale features of the universe 
(see e.g. An essential feature of these inflationary models is their stability, meaning in particular that inflation 
solutions are attractors in the solution space of the relevant cosmological equations (see e.g. 0, @]), at least up until 
inflation ends and the universe becomes radiation dominated. Without this feature, inflation might never have begun, 
i— 1 1 and certainly would not have lasted long enough to affect the large scale structure of the universe. 

Warm inflation is an alternative inflationary scenario in which a small but significant amount of radiation survives 
during the inflationary era due to continuous particle production @, S Q • The coupling between radiation and the 
inflaton field leads to thermal dissipation and fluctuations in the time evolution of the inflaton field. The stab ility 
of the inflationary solutions in warm inflationary models has only been addressed in a limited form previously [1 01 ] . 
and we will present the full stability analysis here. We shall examine conditions under which warm inflation is an 
, ^ , ■ attractor and give conditions for a prolonged period of warm inflation. 

We shall show that the stability of warm inflation can be related to conditions on two parameters describing the 
temperature dependence of terms in the inflaton equation of motion which where not taken into account in the earlier 
stability analysis [Io| . The first condition says that if the dissipation term in the equation of motion falls off too rapidly 
at low temperature then the temperature is driven to zero, and we fall into the conventional inflationarty scenario. The 
second condition limits the temperature dependence of the inflaton potential. In many models large dissipation implies 
large thermal corrections to the potential which prevent warm inflation taking place. This argument is essentially the 
one first presented in Ref. fill ]. Nowadays, we know that there are models where thermal corrections to the inflaton 
potential are suppressed by supersymmetry, and these models may allow warm inflation as an attractor [H, El & 
The duration of the period of inflation is related to a set of slow-roll paramaters which where introduced in [ToT ]. 
We shall re-derive the slow-roll conditions for warm inflation as part of the stability analysis. A well understood 
feature of the slow-roll conditons for warm inflation is that they can be less restrictive than the slow-roll condtions 
. ^ ' for conventional inflation. 

We stress that we are concerned here with the self-consistency of the warm inflationary scenario for given equations 
of motion. We shall not address how the equation of motion for the inflaton field is obtained from non-equilibrium 
thermal field theory. A discussion of the derivation of the equations of motion can be found in a recent review [l2| . 
However, we would like to point out that some of the critisims of warm inflation have been based on models which 
do not satisfy the fundamental stability conditions derived here, and are therefore not inconsistent with the validity 
of warm inflation in general [T(| . 

The stability of warm inflation has consequences for the origin and evolution of cosmological density fluctuations 
[l5l |. In warm inflation, density fluctuations originate from thermal fluctuations [HE?]]- I* 1 particular, the fact that 
the inflationary solution depends on only one parameter means that only one perturbation mode, the curvature 
perturbation, survives on super-horizon scales despite the fact that there are entropy perturbations present on sub- 
horizon scales. We shall give formulae for the spectral indices of the scalar and tensor modes and say a little about 
the tensor/scalar ratio. 
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II. BASIC EQUATIONS 

We start with a flat, homogeneous universe with expansion rate H. The matter content consists of a homogeneous 
inflaton field <p and thermal radiation of temperature T . We restrict attention to the warm inflationary regime where 
T > H, with the radiation close to thermal equilibrium. The evolution of the inflaton is governed by a potential 
V(4>,T) and a damping coefficient T(<p,T), such that the inflaton field satisfies the basic equation 

4> + (3H + r)j> + v <p = o, (l) 

where a comma after a function denotes a derivative. The expansion rate is related to the energy density p by the 
Friedman equation 

3H 2 = 8nGp, (2) 

where G is Newton's constant. 

It is important to realise that the potential appearing in the inflaton equation is the free energy density, rather than 
the potential energy density. The potential energy density is given by the thermodynamic relation V + Ts, where s 
is the entropy density 

s = -V T - (3) 
The total energy density, including the inflaton's kinetic energy density, is therefore 

9 = \i> 2 + V + Ts. (4) 

This includes the contributions of the scalar field and the radiation. It is not always possible to separate the scalar 
and radiation components of the energy density in an unambiguous way. 

The final equation is the one which governs the transfer of energy from the inflaton to the radiation field. This can 
easily be derived from the stress-energy tensor T^ v [l5j , 

7j t „ = Tsu^Uv - Vg^v + V^<jN u (j) - ^(V^) 2 ^, (5) 

where is the radiation fluid 4— velocity vector, g^ v is the metric and V;, is the spacetime derivative operator. 
Conservation of the stress-energy gives 

T{s + 3Hs) = T<P 2 . (6) 

This version of the second law of thermodynamics shows clearly how the friction converts the inflaton's energy into 
heat. We now have a complete set of evolution equations which can be solved for <f> and s given a set of initial 
conditions. 

Inflation is associated with a slow-roll approximation which consists of dropping the leading derivative term in each 
equation. The slow-roll equations are therefore 

Ts = Qtf (8) 
H> - *f (9) 
where the strength of the dissipation is quantified by the parameter Q, 

Q = ^. do) 

Note that, since Eq. (J7JI is first order in time derivatives, any solution to the slow-roll equations has just one constant 
of integration. 

The validity of the slow-roll approximation will depend on the size of a set of slow-roll parameters [l5| . We use the 
following set of 'small' parameters: 

1 fV.,\ 2 _ J_V^£^ 
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An additional pair of parameters describe the temperature dependence, 



V, 



(12) 



We use b in place of the parameter 5 defined in ref [l5[ . 

The parameter b is very important in the theory of warm inflation, ft measures the size of contributions to the 
potential from thermal quantum fields. In ordinary circumstances, we would expect b to be order 1. Consider, for 
example, an inflaton of mass coupled to a set of light scalar fields and temperature T > H > m^. The effective 
potential, 

V (^ T ) = ~^9*T 4 --^mlT 2 +v^), (13) 

where v(<j>) is the effective potential at T = 0. The parameter b for the potential (TIB")) is approximately 2. Values of 
b can be much smaller in supersymmetric theories, where there is a cancellation of leading order thermal corrections 
when the temperature T < As, where As is the supersymmetry breaking scale [181 ] . We shall find limits on b, and 
conclude that warm inflation only takes place when there is a mechanism, like supersymmetry, which reduces the size 
of the thermal corrections to the potential [l!| . 

III. STABILITY ANALYSIS 

We shall consider the consistency of the slow-roll approximation by performing a linear stability analysis to deter- 
mine the conditions which are sufficient for the system to remain close to the slow-roll solution for many Hubble times. 
It may be worthwhile considering first what happens in the alternative cold inflationary scenario (see e.g. Q). The 
slow-roll equation in this case is first order in time derivatives, and the general solution has the form cj> = f(t — to), 
where to is an arbitrary constant. There always exists a homogeneous perturbation which is equivalent to changing 
the value of to, and is hugely important for the existence of density perturbations. Another mode decays on the 
Hubble timescale. It proves convenient to exclude the time-translation mode by using the value of the field as the 
time coordinate, which is possible if the inflaton field has non-vanishing time derivative. We shall follow the same 
procedure for the analysis or warm inflation. 

With the inflaton as independent variable we can rewrite the system of equations in first order form, 



Prime denotes derivatives with respect to <j> and 



where u = <j>. Eqs. (JTJ) and © become 



it 



F(x). (14) 



SH-T-Vjv,- 1 , (16) 



s' = -ZHsu- 1 +T- 1 Tu, (17) 



with the temperature determined implicitly by Eq. (|3|) and H given by Eq. (|2|). We take a background x which 
satisfies the slow-roll equations ([7][5]) and then the linearised perturbations satisfy 

5x' = M(x)Sx - x'. (18) 

where M is the matrix of first derivatives of F evaluated at the slow-roll solution. 

Consider stability first of all. We can express all the components of the matrix M in terms of the slow-roll parameters 
(fT2l . For example, 

cT cT H 

T. s =T, T T s = -T s = -- = cQ-. (19) 
1 6 s s 

where we have used T s = (s ; t) = T/3s. Following a similar procedure for all of the first derivatives gives a final 
expression for M, 

M =(c d) ( 2 °) 
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where 



A = f {- 3{1 + Q) -VTw) (21) 

B = Ij-cQ-^^e+Cl + gjfc} (22) 
C= 5(6-7^1 (23) 



(i + Q) 

(25) 



We require the determinant. 



u 

and the trace, 



d<-1 1/ =^{l2(l + Q) + 3(Q-l)c-6(l + Q)&+ f 3Q ' + 9 ^ 2 + 4 - e + 7-7-77 ) (->(i! 



(i + Q) 2 i + Q; 1 



trM = § {0-4-3(1 + 0)-^}. (27) 

In the cold inflationary case, where Q = b = c = s = 0, the decaying modes have the approximate form Su oc exp(— 3N) 
and 8s oc exp(— AN), where N is the number of e-folds of inflation 

Sufficient conditions for stability of the warm inflationary solution are that the matrix M varies slowly and 

|c|<4-26, 6>0. (28) 

Evidently, the trace is negative and the determinant is positive if these conditions hold. The linear equation therefore 
has two negative eigenvalues and the both eigenmodes decay. The slow variation of M, which allows us to diagonalise 
the linear system, follows if the forcing term in Eq. (|18j) is small, and so we turn to this term next. 

The forcing term in Eq. (|18[) depends on x! . This term is present because the background chosen is not an exact 
solution to the full set of equations. The slow-roll approximation can only be valid when x' is small. If we work with 
time derivatives, the magnitude of x' depends on the the dimensionless quantities u/(Hu) and s/(Hs), and small 
values represent slow variation on the timescale of the Hubble time. 

We only quote the leading terms in e. Starting from the slow-roll Eq. ([9|) and taking the time derivative gives 



H 1 



H 2 1 + Q ' 

By combining the other slow -roll equations (0 and ([8]) we eventually arive at 



(29) 



ii 


1 


Hu 


A 


s 


1 


~Hs~ 


A 



*a+o>»-'fi£^+<«-«to+ I fV>} < 3 ») 

_3(c + g -M - cXl + Q) b + ^ e _ 6|? + 3W-1) J (31) 

where A = 4(1 + Q) + (Q — l)c. The slow-roll approximation requires u << Hu and s << Hs, and sufficient conditions 
for this are 

e, \0\ , \tj\ « 1 + Q; < b « T -^- ! |c| < 4. (32) 

The conditions on e and 77 agree with a previous stabilty analysis [ic|. These are weaker than the corresponding 
conditions for cold inflation, and this fact is a well know feature of warm inflation. 

The physical interpretation of the condition c < 4 is evident from Eq. that radiation must be produced at a 
rate (V oc T c ) exceeding the rate at which radiation is removed by the expansion of the universe (Ts oc T 4 ). The 
condition on b can only be met if there is a mechanism for suppressing thermal corrections to the potential because, 
as we mentioned at the end of Sect [III high temperature thermal corrections would otherwise lead to b « 2. Models 
which include a mechanism for suppressing thermal corrections can be found, for example, in [l3l. [13. [l9| . 
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IV. DENSITY FLUCTUATIONS 



The results which have been obtained as part of the stability analysis are also helpful for analysing various features 
of the density fluctuation spectrum. We therefore take the opportunity, whilst the results are to hand, of giving some 
formulae which might be useful for observational tests of warm inflation. 

The origin of density fluctuations in warm inflationary scenarios is due to thermal fluctuations in the radiation. 
These are coupled to the inflaton as a consequence of the friction term in the inflaton equation of motion, and their 
amplitude is fixed by a fluctuation-dissipation theorem. This means that both entropy and curvature perturbations 
must be present. However, on length scales larger than the horizon, we know from the stability argument that 
the coupled inflaton plus radiation system approaches the slow-roll solution which has only one free parameter. 
Consequently, on large scales only the pure curvature perturbation survives. This has been confirmed in particular 
models by solving the full set of density fluctuation equations numerically [Til ] . 

Even though the entropy perturbations decay on large scales, they can sometimes leave behind an impression on the 
curvature fluctuations. If the friction term depends on temperature, then the entropy and curvature fluctuations on 
sub- horizon scales become coupled. The situation is similar to the sympathetic oscillations of a double pendulum [20| . 
When the curvature fluctuations 'freeze-out', the amplitude of the sympathetic oscillation may be anywhere between 
zero and its maximum value, leading to oscillations in the wave-number dependence of the curvature fluctuation 
spectrum. The amplitude given below therefore has only limited use when b and c are non-zero and refers only to the 
envelope of these oscillations. 

The thermal fluctuations produce a power spectrum of scalar density fluctuations of the form [TBI, [2l[ , 



The spectral index n s is defined by 



-Ps = ^^f(l + Q) 1/2 - (33) 



1 dlnVs CiA\ 

ns ~ l = - aEfc-' (34) 



evaluated when the amplitudes 'freeze out'. To leading order in the slow-roll parameters we can take the freeze-out 
time to be the horizon crossing time when k — aH, and then 

We can use eqs. (I29H31[) to obtain 

1 1 f 3(2Q + 2 + 5Qc)(l + Q) L 9Q + 17-5c 9Q + l a 3Qc-6-6Q + 2c 1 
ns ~ 1 = A 1 Q b - T+Q 6 TTqP 1 + Q V j (36) 

If we consider b = c = 0, then important limits include the strong regime of warm inflation, Q >> 1, 

" s - 1 = -^-^" + ^"- < 37 > 

This agrees with a partial result in [22| and the full result in (l5| . In the weak regime of warm inflation, Q « 1, 
thermal fluctuations lead to the spectral index 

ns-l = ~e-- A P+^n. (38) 

Previous results for the weak regime, though expressed in a less useful form, can be found in Refs (23l. [24|. Finally, 
the case Q << 1 and c = 3 is important because it corresponds to a class of warm inflationary models where the 
friction coefficient T has been calculated [25| . 

n s -l = -2e- P. (39) 
The tensor modes have the same amplitude as they do in the cold inflationary models, (2|| 



T T = 8nGH 2 . 



(40) 



6 



The spectrum for the tensor modes is simply 



n T -l = - * e. (41) 



Unlike in the cold inflationary scenario, the tensor-scalar amplitude ratio cannot be expressed in terms of slow-roll 
parameters. Instead [2^ |. 

Vt_ 2e H 

V S {1 + QfT- K ' 

Since T > H for warm inflation, the tensor-scalar ratio is likely to be smaller than 1 — tit- Tensor modes are strongly 
suppressed relative to the scalar modes on the strong regime of warm inflation Q » 1, but they could be significant 
in the weak regime of warm inflation. 

We conclude this section by finding the lower limit on the friction term which is required for warm inflation. We 
shall express this limit in terms of Q = T/3H. Re- write the scalar amplitude Eq. (|33|) as 

Vs^%^{l + Q) 1/2 (43) 
u T 

The first factor can be replaced using the slow-roll equation Eq. ([8]) and the potential (|13p . and we obtain 

r-./ « v'vw^V". ,44, 



The condition for warm inflation T > H is therefore 

Q > 9*Vs- (45) 

Cosmic microwave background observations give a scalar power spectrum of 1 x 10 -10 on large scales, therefore very 
small amounts of dissipation can result in warm inflation. 



V. CONCLUSION 



We shall recapitulate the main points of this paper. There are conditions on six of the parameters defined in Sect. 
[11] for the possibility of a stable period of warm inflation in the early universe: 

• The parameter Q which measures the strength of the friction term must satisfy 

Q > 9*V S . (46) 
where is the effective particle number and Vs is the scalar perturbation power spectrum on large scales. 

• The parameters which describe the inflaton dependence of the effective potential and friction term satisfy 

e«l + Q, M«l + Q, \/3\«l + Q. (47) 



• The temperature dependence of the potential and the friction term is restricted by 

\b\«Y^Q, N<4. (48) 

The condition on b implies that warm inflation is only possible when a mechanism, such as supersymmetry, 
reduces the size of the thermal corrections to the potential. 

Models of elementary particles exist in which these conditions can be satisfied. The most convincing of these models 
use a combination of supersymmetry and a two-stage decay process, where there are no direct coupling between the 
inflaton and the radiation and all thermal effects are supressed by factors of T/Ag, where Ag is the supersymmtry 
breaking scale [H Q [lj] . 

When the conditions listed above are satisfied, then the solutions to the equations of motion approach a slow-roll 
approximation during inflation. As a result, large scale density perturbations have only one degree of freedom, which 
we identify as the curvature perturbation. (Entropy perturbations can only be introduced by adding extra degrees of 
freedom to the system.) 
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